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Using a weak-disorder scheme and real-space renormalization-group techniques, we obtain analytical results 
for the critical behavior of various q-state Potts models with correlated disordered exchange interactions along di 
of d spatial dimensions on hierarchical (Migdal-Kadanoff) lattices. Our results indicate qualitative differences 
between the cases d — d\ — 1 (for which we find nonphysical random fixed points, suggesting the existence 
of nonperturbatfve fixed distributions) and d — di > 1 (for which we do find acceptable perturbartive random 
fixed points), in agreement with previous numerical calculations by Andelman and Aharony. We also rederive a 
criterion for relevance of correlated disorder, which generalizes the usual Harris criterion. 
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I. INTRODUCTION 

The effects of disorder on the critical properties of statis- 
tical models have been the subject of much work in the last 
decades. In the context of random interactions, Harris [|lj de- 
rived a heuristic criterion to gauge the relevance of uncorre- 
cted disorder to the critical behavior, which is predicted to re- 
main unchanged if the specific-heat exponent a of the under- 
lying pure system is negative. If a > 0, disorder becomes rel- 
evant and, in the language of the renormalization group (RG), 
one expects a flow to a new fixed point (characterized by a 
nonzero- width fixed distribution of the random variables). 

It later became clear that the Harris criterion must be gener- 
alized in a number of situations [|[ ||, |],|^, §], since a is not al- 
ways identifiable with 0, the crossover exponent of the width 
of the distribution of the disorder variables. In particular, ran- 
dom variables correlated along di of the d spatial dimensions 
give rise to the scaling relation | 



(1) 



■dif, 



where v is the correlation-length exponent of the pure system. 
Using a real-space RG approach based on numerical calcu- 
lations [R], Andelman and Aharony [|j] investigated various 
g-state Potts models with random exchange constants, finding 
qualitative differences between the cases d — d\ > 1 (which 
yields finite-temperature fixed distributions) and d — d% = 1 
(which embodies the McCoy-Wu model [Q] and yields an 
"infinite-disorder" zero-temperature fixed point). An intuitive 
illustration of the special role of the d — d± = 1 case is that, 
for any infinitesimal concentration of zero bonds (with a suit- 
able assignment of the random interactions), the system would 
break into non-interacting (d— 1) -dimensional structures, and 
the RG flows would be redirected to the pure fixed point of 
the corresponding system in d — 1 dimensions. 
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In the present paper, we use a (perturbative) weak-disorder 
[§, |Io|] real-space RG scheme to analyze the critical behavior 
of g-state Potts models with correlated disordered exchange 
interactions on various hierarchical lattices, whose exact re- 
cursion relations are equivalent to those produced by Migdal- 
Kadanoff approximations for Bravais lattices. Using this 
weak-disorder scheme, we obtain analytical results by trun- 
cating the recursion relations for the moments of the disorder 
distribution (which are supposed to remain sufficiently small 
under the RG iterations). All calculations are performed in 
the vicinity of <f> = 0, in a region where disorder is rele- 
vant. Depending on the diference between the dimensionality 
of the system (d) and the number of dimensions in which dis- 
order is correlated (di), we distinguish two possibilities: (i) 
For d — di = 1, the weak-disorder scheme produces a non- 
physical fixed-point probability distribution, characterized by 
a negative variance, which suggests the existence of a nonper- 
turbative ("infinite-disorder") fixed-point; (ii) For d — d\ > 1, 
the scheme yields a physically acceptable perturbative fixed- 
point distribution. Although obtained by an alternative ap- 
proach, the main results of this paper are in agreement with 
the numerical findings of Andelman and Aharony [Q]. 

The outline of the paper is as follows. We first rederive 
Eq. ([[]), and obtain a criterion for relevance of correlated dis- 
order involving the number of independent random variables 
in the unit cell of the lattice and the first derivative of the re- 
cursion relations at the pure fixed point. This is done in Sec. 
U In Sec. H we consider g-state Potts models on various 
hierarchical lattices with d — d\ — 1. Using a weak-disorder 
scheme, we obtain a new (random) fixed point for q larger than 
a characteristic value qo, where disorder becomes relevant. As 
in a previous publication JlQ[], this fixed point is located in a 
nonphysical region of the parameter space, suggesting that a 
nonperturbative fixed point must be present. In Sec. [v| we 
study a similar problem with di = 1 and d = 3. In this case 
we obtain a physically acceptable, finite-disorder fixed point, 
for q > qo, as in the fully disordered model studied by Derrida 
and Gardner (although in our case the usual Harris crite- 
rion is not satisfied). In Sec. |v|, we consider an Ising model 
(q = 2) on a diamond lattice with b = 2 bonds and I branches 
(where I, instead of q, is the control parameter), which con- 
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stitutes another example of a d — d% = 1 system. As in Sec. 

weak disorder again predicts a nonphysical random fixed 
point. In the final section we give some conclusions. 



II. CRITERION FOR RELEVANCE OF CORRELATED 
DISORDER 

Following Andelman and Aharony [jjj], we consider a 
e?-dimensional bond-disordered model in which the dis- 



order variables are correlated along d\ spatial directions. 
We assume that, under renormalization with a length 
rescaling factor b, the model satisfies a recursion relation 
R (xi, X2, ■ ■ ■ , x n ), connecting n = b d ~ dl independent (and 
identically distributed) random variables to a renormalized 
variable x'. (In this paper, these variables are related to re- 
duced exchange couplings.) Defining the deviations Si = 
Xi — x c , where x c — R(x c , x c , . . . , x c ) is the critical fixed 
point of the pure system, we expand R in a Taylor series about 
x r to write 
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and similarly for the higher powers of e'. Averaging over the random variables we get 
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and corresponding expressions for the higher moments of the 
deviations. Since (e) is a measure of the distance to the fixed 
point, it plays the role of temperature. On the other hand, (e 2 ) 
is a measure of the strength of disorder. 

The critical behavior of the model is related to the eigenval- 
ues of the matrix 



M r , = 



d(e s ) ' 



(6) 



evaluated at the fixed point. It is clear that the set of recursion 
relations for the moments of the deviations always has a pure 
fixed point (e) = (e 2 ') = ■ • • = 0. At that point, it can be 
shown [|fl]] that M rs is a triangular matrix, and that its two 
largest eigenvalues are given by 



Ai 
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and 



A 2 = 



die'- 



die 2 ) 



pure 
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Assuming that, for all i and j, 

dR dR 

dxi _ dxj 



w, 



and invoking the usual scaling hypotheses 

Ai = b yt and A 2 = A? = 



(9) 



(10) 



which define the thermal exponent y t and the crossover expo- 
nent <f>, we get 



(jnjt = 2y t - (d - di) 
Then, using the hyperscaling relation 



a = 2 = 2 

Vt 



d\nb 
hx{nw) ' 



we obtain 



di d-di di 

a H = — a + 2—, 

yt d d 



(ID 



(12) 



(13) 



which clearly shows that the Harris criterion (<t> — a > 0) is 
not satisfied in the presence of correlated disorder. As l/y t 



FIG. 1: (a) The diamond hierarchical lattice (with b — 2 and 1 — 2). 
(b) The necklace hierarchical lattice (with 6 = 2 and 1 = 2). 



is usually identified with the correlation-length exponent v, 
this last result is equivalent to Eq. ([j]). It also shows that, for 
d\ > 0, the crossover expoent is larger than a, which indi- 
cates that correlated disorder induces stronger (geometrical) 
fluctuations than uncorrelated disorder. 

The general criterion for relevance of disorder is <p > 0, 
that is, 



a > -2- 



di 



d— di 

From Eqs. (0)-(§), this is equivalent to 



nw > 1. 



(14) 



(15) 



This last result was also derived in a different context by 
Mukherji and Bhattacharjee [p|] and generalizes a criterion 
pointed out by Derrida et al. [M. 

In the case of the fully disordered system analyzed by Der- 
rida and Gardner [^], for which d\ — 0, the requirement in 
Eq. (|l4|) turns out to be equivalent to the usual form of the 
Harris criterion (a > 0). 



III. POTTS MODELS WITH CORRELATED DISORDER: 

d - di = 1 CASE 

The successive generations of a hierarchical lattice are ob- 
tained by replacing an existing bond in the previous genera- 
tion by a unit cell of new bonds in the next generation. In 
Fig. [jja), we show the first two stages of the construction 
of the simple diamond lattice (with 6 = 2 bonds and I = 2 
branches). The necklace hierarchical lattice, with 6 = 2 bonds 
and / = 2 branches, is illustrated in Fig. |](b). 

We now consider a g-state Potts model, given by the Hamil- 
tonian 



(16) 



where the sum is over nearest-neighbor sites on a hierarchical 
lattice, the spin variables cr, assume q values, S is the Kro- 
necker symbol, and {Jjj > 0} is a set of independent and 
identically distributed random variables. Instead of consid- 
ering a fully disordered arrangement of interactions, we look 
at correlated disorder, either along layers (see Figs. |](a) and 
^(c)) or along branches (see Figs. ||(b) and^(d)) of the hierar- 
chical structure. 

Introducing the more convenient variable x. L — exp(/3Jj), 
where (3 is the inverse absolute temperature, it is straightfor- 
ward to decimate the internal degrees of freedom to obtain 





FIG. 2: Correlated distribution of random interactions on diamond 
and necklace hierarchical lattices. 



(Migdal-Kadanoff) recursion relations. In this section we con- 
sider the following models: 

A. random layered diamond lattice, Fig. ^|(a), whose re- 
cursion relation is 



x' = Ra (xi,x 2 ) 



x\x 2 + q — 1 
%\ + x 2 + q - 2 



(17) 



B. random branched diamond lattice, Fig. ||(b), with re- 
cursion relation 



x 1 = Rb(xi,x 2 ) = 



I + q - 1 \ / xj + q-1 



2xi+q-2J \2x 2 + q-2 



(18) 

C. random layered necklace lattice, Fig. ||(c), with recur- 
sion relation 



X = R C (X!,X 2 ) = 2 1 _^ 2 '2 9 ~ 



x\ + x\ + q — 2 



(19) 



D. random branched necklace lattice, Fig. ^(d), with re- 
cursion relation 



x = R D (xi,x 2 ) = 



x\x 2 + q — 1 
2x\X 2 + q — 2 



(20) 



Notice that in all these models disorder is correlated along 
only one spatial direction (di — 1), while the effective dimen- 
sion is d = 2. According to Eq. (p^j), we then expect disorder 
to be relevant for a > — 2. 

We now write x' = x c + e' and xi — x c + Si, to perform 
Taylor series expansions about the critical point of the uniform 
systems, given by x c = R(x c , x c ). For all of these models, 
with n = 2 independent values of the exchange parameters 
(along either layers or bonds), it is straightforward to write 
the recursion relation 



w (ei +e 2 ) + m (ej + ef) + f (eiej + e\e 2 ) 



pei£ 2 + ce\e\ + k(el + el) + a (ej 



(21) 



where w, m, p, f, c, k and a are model-dependent Taylor coef- 
ficients (that depend on the topology of the particular models 
illustrated in Fig. ^; see Sec. |l|). 
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The weak-disorder approximation [T(]] consists in as- 
suming that 



and in general 



(e) ~ (e 2 ) ~ A, 



<e 3 ) ~ < £ 4 ) ~ A 2 



(22) 



(23) 



(24) 



where (■ ■ •) is a quenched average and A is a suitable small 
parameter. Within this approximation, we can use Eq. ([ll]) to 
write recursion relations for the moments of the deviation, up 
to second order in A, 

(e') = 2w(e)+ P (e) 2 + 2m (e 2 ) + 2/ (e) (e 2 ) 

+ c( £ 2 ) 2 + 2fc( £ 3 ) +2a(e 4 ), (25) 



= 2w 2 (e) 2 + 2w 2 (e 2 ) + 4w(m + p) (e) (e 2 ) 

+ (2m 2 + ifw +p 2 ) (e 2 ) 2 + iwm(e 3 ) 
+ (4wk + 2m 2 ) (e 4 ) , 



and 



= 3w(e) (e 2 )+3(m+p) {e 2 f 
+ w(e 3 ) +3m(e i ) , 

(s' 4 ) = 3w 2 (e 2 ) 2 + w 2 (e 4 ) . 



(26) 



(27) 



(28) 



It is easy to see that there is always a non-random fixed 
point, 



(e) = (e 2 ) = (e 3 ) = (e 4 ) = 0, 



(29) 



associated with the critical behavior of the pure model. As we 
pointed out in the previous section, this fixed point becomes 
unstable with respect to disorder for 2w 2 > 1. This can also 
be seen by an inspection of the asymptotic behavior of Eq. 
([26]), which shows that, up to order A, the renormalized sec- 
ond moment depends only on (e 2 ), with the coefficient 2w 2 . 
Thus, we expect the onset of a random fixed point at a critical 
value go of the number of Potts states. From the expression 



x c R(x c ,Xc) 



(30) 



for the pure fixed point, we can express q as a function of x c 
and, using the condition 2w 2 = 1, determine the critical value 
x c (on). For both diamond structures displayed in Figs, 0^a) 
and 0(b), we have 



q = (y/xc- !) ( x c - 1) : 



(31) 



and x c (qo) — 2.15127. . ., which leads to qo = 0.53732 . . 
For both necklace structures in Figs. ||(c) and||(d), we have 

q=(x c - 1) (x 2 - 1) , 



coefficient model A model B model C model D 

a -0.00926 0.00917 -0.92623 0.02894 

e 0.08549 0.00016 1.38173 0.07163 

k 0.04676 -0.01302 0.25648 -0.02801 

/ -0.05370 0.00608 -0.33156 -0.04706 

p 0.65117 0.23242 1.56929 0.53634 

TABLE I: Coefficients of the weak-disorder expansion for the mod- 
els in Fig. |i| 



moment model A model B model C model D 
(e) /Aq -14.904 1.0208 -4.4401 0.34798 



/Aq 



-16.170 -11.434 -1.8791 -2.6575 



/(A 9 ) 2 



1444.5 325.73 46.390 39.946 



(e 4 ) /(Agf 784 



41 392.21 10.593 21.187 



TABLE II: Moments of the deviations defining the random fixed 
points of the models in Fig. ^ according to the weak-disorder ex- 
pansion. 



with x c (qo) = 1.46672..., which also leads to q = 
0.53732 . . .. Disorder is predicted to be relevant for q > qo. 
We now introduce the small parameter 



A = x c (q)-x c 0? o ) 



dx c 
dq 



(q-qo) 



dx c 
dq 



Ag, (33) 



to investigate a q-state Potts model in the immediate vicinity 
of the characteristic value qo. It should be pointed out that, 
as the symmetry of the order parameter is one of the factors 
expected to determine the universality class of the models, Aq 
is the appropriate parameter to consider. However, A is more 
convenient for the algebraic manipulations. From inspection 
of Eqs. (^5|)-(^), we see that, up to first-order terms in A, 
coefficients w and m are written as 



w 



wiX 



and 



mo + TOiA. (34) 



It is straightforward to calculate w\ = 0.13325 . . ., for the 
diamond structures, and w-y — 0.39088 . . ., for the neck- 
lace structures. Also, we have mo — —0.19088... and 
mi = 0.19865..., for model A; m = 0.01849... and 
mi = 0.00758 . . ., for model B; m = -0.48935 ... and 
mi = 1.22433 . . ., for model C; and m = 0.02711 ... and 
mi = 0.02027 . . ., for model D. In order to obtain the remain- 
ing coefficients, it is enough to keep the zeroth order term in 
A (see the values, up to five digits, in Table |). 

We are finally prepared to obtain, up to lowest order in Aq, 
the non-zero values of the moments at the random fixed point. 
By substituting the weak-disorder assumptions, Eqs. ( |22| ) and 
(|23|), into Eqs. (^)-(^8j), and then imposing consistency be- 
tween equal powers of Aq, we obtain the leading terms for 
fixed values of the moments as listed in Table |l[ 

In order to perform a linear stability analysis about the fixed 
(32) points, we have to calculate the eigenvalues Ai to A4 of the 
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matrix 



M r 



d(e s ) 



As it should be anticipated from universality, it turns out that 
the eigenvalues (and so the critical exponents) are the same 
for all models A to D. We always have two eigenvalues, A3 
and A4, whose absolute values are smaller than unity. About 
the pure fixed point, we have 



A^ p) = \/2 + 0.31018A g , 
0.43866Ag, 



A« = 1 



(35) 
(36) 



with a specific heat exponent 

a p = -2 + 2.5314lA g . 
At the random fixed point we have 



= V2 + 0.83670Ag, 



A 2 r) = l-0.43866A g , 



which lead to the exponent 



a t = -2 + 6.82843Ag. 



(37) 
(38) 

(39) 



From Eq. (|3q), we see that disorder becomes relevant for 
Aq > 0. Thus, as shown in Table||, the weak-disorder expan- 
sion gives negative (and thus nonphysical) values of the sec- 
ond moment at the random fixed point for all models A to D. 
This suggests that the random fixed point in these systems (for 
which d — di = 1) is nonperturbative, in agreement with nu- 
merical calculations [Q] that predict an infinite-disorder fixed 
point. Another odd feature of the weak-disorder results is that 
the predicted value of the specific-heat exponent in the pres- 
ence of disorder (a r ) is larger than the corresponding quantity 
(a p ) for the pure model, in disagreement with the general be- 
lief that disorder should weaken the transition. 
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IV. A POTTS MODEL WITH CORRELATED DISORDER: 

d - di > 1 CASE 



In order to examine the d — d\ > 1 case, we now consider 
a Potts model on a necklace hierarchical lattice shown in 
Fig. ||, with d = 3 and di = 1. The unit cell contains n = 4 




FIG. 3: The hierarchical lattice with d = 3 and d\ = 1 considered 
in Sec. [iv] 

independent random variables and, in terms of the variables 
X{ = exp(/3 Ji), the recursion relation is given by 



R(Xl,X 2 ,X 3 ,X i ) = ; ; -. (40) 



X\X2 + X3X4 + q - 2 



Following the same steps as in Sec. 5BL we have 



q=(x c - 1) {x\ - 1) , 



(41) 



go = 4 + 2\/2, and x c (qo) = 1 + \/2. Performing again 
the weak-disorder expansion (and truncation), and taking the 
average over the disorder variables, we obtain the set of recur- 
sion relations 



<e') = Aw (e) + 2(pi + 2 P2 ) (ef + Am (e 2 ) + 4(/i + 2/ 2 ) (e) (e 2 ) 
+ 2( Cl + 2c 2 )(e 2 ) 2 + 4fc(e 3 )+4a(e 4 ), 



(42) 



and 



./ 3 



= 12w 2 {e) 2 + 4w 2 (e 2 ) + 8w(3m + Pl + 2p 2 ) (e) (e 2 ) 

+ [12m 2 + 8w(h + 2/ 2 ) + 2(p 2 + 2p 2 )} (e 2 ) 2 
+ 8wm(e 3 ) + (8wfc + 4m 2 ) (e 4 ) , 



\ = 9w (e) (e 2 ) + 3(3m + Pl + 2p 2 ) (e 2 f + w (e 3 ) + 3m (e 4 ) 



(e' 4 ) = 9 W 2 (e 2 ) 2 + ^ 2 ( £ 4 ) 



(43) 
(44) 

(45) 
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It should be noted that, due to the smaller symmetry of the 
lattice, we now have a larger set of coefficients. Also, notice 
that in this case go is determined from the condition Aw 2 = 1. 
About the critical value g , and to leading order in Aq, we 
have 



w 



1 17V2-24 



Aq 



and 



\f2-2 133- 94^ 
8 + 16 



Aq. 



(46) 



(47) 



The values for the remaining coefficients are listed in Table 


The moments of the deviations at the random fixed point 
are written as 



<^ 2 > 
<^ 4 > 



i (5 - 3^) Aq, 
I (a - V2) Ag, 

(95V2- 128^ (Aq) 2 
9 - Ay/2] (Aq) 2 . 



7 

3 , 

49 
6 

49 



(48) 



Performing a linear stability analysis about the pure fixed 
point we obtain 



A< p) 



2+ (m/2-24) Aq, 



A 



(p) 



1+ (17V2-24J Aq, 



with a specific-heat exponent 

3 17V2- 24 . 
a ^- l+ 2 In 2 Agl 

while about the random fixed point we have 



A 



(r) 



2 - - (92 - 65\/2) Ag, 



A^ r) = 1 - (l7\f2- 24) Ag, 



with 



Or 



-1 



3 92 - 65V2 
14 m~2 



Ag. 



(49) 
(50) 

(51) 

(52) 
(53) 

(54) 



These results show that once more disorder becomes rele- 
vant for Ag > 0, but now we obtain a positive (and thus physi- 
cally acceptable) value of the second moment of the deviations 
at the random fixed point. We also have a r < a p . So, as in the 
fully disordered model (d± = 0) studied by Derrida and Gard- 
ner [^J, and in agreement with numerical calculations [[5J, the 
weak-disorder scheme predicts a (perturbative) finite-disorder 
fixed point, with values of the critical exponents continuously 
approaching those of the pure model as Ag — > 0. 




I branches 

FIG. 4: A diamond hierarchical lattice with b 
branches. 



2 bonds and I 



V. AN ISING MODEL WITH CORRELATED DISORDER 

The set of recursion relations given by equations ( p5| ) to 
(|28|), with a suitable redefinition of parameters, can also be 
used to analyze an Ising model on a more general diamond 
structure with 6 = 2 bonds and I branches, and correlated dis- 
ordered ferromagnetic exchange interactions along the layers 
(see Fig. Q). For this structure we also have d — di = I. While 
in the Potts models we have a natural parameter, q, for varying 
a, we now change the topology of the lattice, by varying I, to 
obtain the same effect. 

Using the standard Ising Hamiltonian, 



(i,3) 



(55) 



with <Tj = ±1, and introducing the more convenient transmis- 
sivity variable ti = tanh j3 Ji, the decimation of the interme- 
diate spins leads to the recursion relation 

t' = R((h,t 2 ) = tanh{/tanh _1 (tit 2 )}. (56) 

As in Sec. Ill, we now write t' = t c + e and ij = t c + £.- L , 
where 



t c — Ri(t c , tc) 



(57) 



is the critical transmissivity of the uniform model. We then 
perform quenched averages, and use the weak-disorder as- 
sumption, to obtain Eqs. (^5j) to (p8|). 

The critical parameters for relevance of disorder, Iq — 
1.44976 ... and t c (l Q ) = 0.79951 . . ., come from Eqs. @ 
and (|l5|). The small parameter A can be chosen as 



A = t e (I) - U (Jo) = -Jr 
dl 



AJ. (58) 



Again we use A as a convenient parameter for algebraic ma- 
nipulations, although Al is the physically relevant variable. 
The Taylor coefficients in Eqs. ( p5| ) to ( |28| ) are given by 
w = V2/2 - 0.54522A, m = -0.49698 - 0.65422A, a = 
0.11520, c = -1.64903, k = -0.12543, / = -1.61924, and 
p = —0.10953. We then calculate the leading values of the 
moments at the random fixed point, 

(e) = -0.64971AZ, 
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pi 


P2 


Cl 


C2 


/l 


/2 


ft 


a 




1^-1 


109\/2-144 


27^2-38 


25-18^2 


11-8V2 




7^2-10 


4 


X 2 


32 


32 


lfi 


lfi 


lfi 


fi4 



TABLE III: Values of the weak-disorder coefficients for the model in Sec. UV 



<^ 2 > - 
<^ 3 > = 
<^ 4 > = 



-0.27076AZ, 

-0.30084(AZ) 2 , 

+0.21993(AZ) 2 . 



(59) 



A linear stability analysis leads to the eigenvalues A^ p -* = 
\/2+0.71884AZand A^ p) = 1+ 1.01659 AZ, for the pure fixed 
point, and A^ = V2+1.20537AZ and A^ r) = 1-1.01659AZ, 
for the random fixed point. From these values, we see that dis- 
order is relevant for AZ > 0, but we again have ^e 2 ) < in 
this case. 

We then obtain the specific heat critical exponents 



a p = -1.07163 



2.51471AZ 



and 



-1.07163 + 5.56379 AZ. 



(60) 



(61) 



For AZ < 0, which corresponds to a < —1.07163 . . ., the 
pure fixed point is stable and the random model displays the 
same critical behavior as its pure counterpart. For AZ > 0, 
which corresponds to a > —1.0713 . . . (yielding again a T > 
a p ), we anticipate a novel class of (random) critical behavior, 
but the fixed point must be nonperturbative, as suggested by 
the nonphysical character of the weak-disorder results. 



VI. CONCLUSIONS 

We have used a weak-disorder scheme and real-space 
renormalization-group techniques to look at the effects of cor- 
related disorder on the critical behavior of some (/-state Potts 
models with correlated disordered ferromagnetic interactions 
along eZi out of d spatial dimensions. We have written ex- 
act recursion relations on diamond and necklace hierarchical 
structures, which are equivalent to Migdal-Kadanoff approxi- 
mations for the corresponding Bravais lattices. 

The weak-disorder scheme leads to analytical results by 
truncating the recursion relations for the moments of the dis- 
tribution function. We first used scaling arguments to rederive 
a general expression for the Harris criterion to gauge the rel- 
evance of disorder (and show that it is related to the number 
of independent random variables in the unit cell of the lattice 
and the first derivative of the recursion relations at the pure 
fixed point). We then performed a number of calculations to 
compare with numerical findings by Andelman and Aharony. 



For q-state Potts models on various hierarchical lattices 
with ferromagnetic random exchange interactions correlated 
along di = 1 out of d = 2 dimensions, we obtained a new 
(random) fixed point for q larger than a characteristic value qo, 
where disorder becomes relevant. This fixed point, however, 
is located in a nonphysical region of parameter space, which 
suggests that a nonperturbative (infinite-disorder) fixed point 
must be present (as pointed out by the calculations of Andel- 
man and Aharony). For a g-state Potts model on a diamond 
lattice with di = 1 and d — 3, we obtained a physically ac- 
ceptable, finite-disorder fixed point, for q > qo, as in the fully 
disordered model analyzed by Derrida and Gardner (although 
in our case the usual expression of the Harris criterion is not 
fulfilled). Also, we considered an Ising model (q = 2) on a 
diamond lattice with 6 = 2 bonds and Z branches (where Z, in- 
stead of q, is the control parameter), which is another example 
of a d — d\ = 1 system. Again, the weak-disorder expansion 
predicts a nonphysical random fixed point. 

To summarize the results of this paper, we point out that, in 
the vicinity of the point where disorder becomes relevant, the 
weak-disorder scheme always produces a perturbative random 
fixed point, but there are two distinct possibilities, depending 
on the difference between d and d\. (i) If d — d\ = 1, the per- 
turbative fixed point is characterized by a negative variance, 
and is thus nonphysical, suggesting the existence of another, 
nonperturbative fixed point; (ii) If d — d\ > 1, the scheme 
predicts a physically acceptable perturbative fixed point. It 
should be mentioned that this same picture holds for fairly 
general hierarchical lattices, in particular those with noniterat- 
ing bonds, as considered by Griffiths and Kauffman [ 12 1. Fur- 
thermore, in the case of the quantum Ising model with bond 
disorder, which corresponds to the extreme-anisotropy limit of 
the two-dimensional McCoy-Wu model (d — d\ = 1), Fisher 
Jl3| ] was able to obtain a (presumably exact) fixed-point prob- 
ability distribution with infinite variance. It is certainly inter- 
esting to investigate whether similar conclusions still hold for 
other models (as the problem of directed polymers in random 
environments [|]|) on either hierarchical or Bravais lattices. 
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